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Abstract 

In this note, we study constraints on primordial non-Gaussianity in DBI Galileon models in which an 
induced gravity term is added to the Dirac-Born-Infcld (DBI) action. In this model, the non-Gaussianity 
of orthogonal shape can be generated, which is slightly favoured by WMAP nine-year (WMAP9) data. 
We provide a relation between theoretical parameters and orthogonal/equilateral non-linear parameters 
using the Fisher matrix approach. In doing so, we include the effect of the CMB transfer functions and 
Planck noise properties by employing the recently developed SONG code. The relation is also shown in 
the effective theory language so that it can be applied to general single-field models. Using the bispectrum 
Fisher matrix for Planck, we derive forecasts for constraints on the theoretical parameters. We find that 
Planck can rule out the conventional DBI inflation model at a statistically significant level, provided that 
it measures the same central values as WMAP9 for equilateral and orthogonal non-Gaussianities. 

1 Introduction 

Primordial non-Gaussianity of the curvature perturbation provides valuable information on the physics in the 
very early Universe pQ. Non-linearity of quantum fluctuations during inflation gives rise to a bispectrum in 
the Cosmic Microwave Background (CMB) temperature anisotropies that peaks for equal-side triangles [2J. 
The most popular model for the equilateral type non-Gaussianity is the Dirac-Born-Infeld (DBI) inflation 
model [3j. In the DBI inflation model, the inflaton is identified as the position of a D3-brane in a higher 
dimensional spacetime. The DBI action that describes the motion of the brane is a non-linear function of 
the kinetic term of the inflaton, which leads to the non-Gaussianity of quantum fluctuations. While the 
original DBI inflation model considered only the motion and fluctuations in the radial direction, we can 
consistently take into account the dynamics and fluctuations in the angular directions [1] (see also [6] ) . 
Although the original single field DBI inflation model [3] is under strain from an additional requirement 
related with the compactification scheme of the string theory [8] , this can be evaded in multi- field DBI 
inflation models and the shape of the bispectrum remains as in the single-field model [1]. 

Recently, a natural extension of the DBI inflation model has been obtained by adding an induced 
gravity term [1U\ 111) . This leads to the quartic Galileon Lagrangian [T2J when the motion of the brane 



is non-relativistic. Thus this model is dubbed as the DBI Galileon model [13]. This is one of a very few 
models where the non-Gaussianity of orthogonal shape can be generated [lUj . The orthogonal shape of 
non-Gaussianity was originally discovered in the context of the effective theory of inflation [T3], which has a 
minimum overlap between local and equilateral non-Gaussianities [TS]. In the WMAP nine- year (WMAP9) 
data, a hint was found that the orthogonal type non-Gaussianity could be non-zero at 2a level when the 
equilateral non-Gaussianity is included in the parameter space [9]. This result is far from conclusive and it 
needs to be confirmed by the Planck satellite. 

In this note, we provide a relation between the equilateral and orthogonal templates, parametrised 
respectively by the non-linear parameters and f^^ 1 , an< ^ the theoretical parameters in the DBI Galileon 
models by properly taking into account the CMB transfer functions and Planck noise properties. We also 
derive the constraints on /^li/nl^ ex P ec t e d from Planck, including the possibility of the simultaneous 
presence of the two shapes in the data. We then provide forecasts for the constraints on the parameters in 
the DBI Galileon model. We also present these constraints in the effective theory language so that they can 
be easily applied to more general single-field models. 

This paper is organised as follows. In section 2, we summarise the prediction of non-Gaussianity in the 
DBI Galileon models. In section 3, we present the equilateral and orthogonal templates and discuss the 
overlap between theoretical bispectrum shapes and these templates. In section 4, we apply the WMAP9 
results to obtain constraints on the theoretical parameters. We study forecasts for Planck in section 5. We 
compute the bispectrum Fisher matrix for Planck and use it to provide a relation between observational 
templates and theoretical parameters. Similarly, we provide Planck forecasts for the two parameters in the 
DBI Galileon model and those in the effective theory. Section 6 is devoted to the conclusion. 



2 Non-Gaussianity in DBI Galileon model 

In this section, we summarise the set-up of the DBI Galileon model and its predictions for the non- 
Gaussianity by following Refs. [101111]. We consider a D3-brane with tension T3 evolving in a 10-dimensional 
geometry described by the metric 

ds 2 = h-V^y") g^dx^dx" + ^/ 2 (y*) G u (y K ) dy 1 dy J = H AB dY A dY B , (1) 

with coordinates Y A = {x^, 7/}, where \i = 0, . . . 3 and I = 1, ... ,6. The induced metric on the 3-brane is 
given by 

1 ^ = H AB d^Y [ i ) d u Y { l ) , (2) 

where the brane embedding is defined by the functions Y^Mx* 4 ), with x^ 1 being the spacetime coordinates on 
the brane. We choose the brane embedding as Y,A = (x^, ip 1 (x 11 )). Then the induced metric can be written 



as 

7m- = h ~ 1/2 (<V + h Gijd fl <p I d v <p J ) . (3) 
The action in the DBI Galileon model is given by 



s = hrx 



Ml , „ r , M 2 



2 v J LJJ 2 
where Mp and M are constant mass scales and 



^g-R[g] + ^-^Rfr] + 



brane 



(4) 



£branc = (VV - l) - V{<?) . (5) 

The second term in the action is the induced gravity term, which is absent in the conventional DBI inflation 
model. Here, we have introduced the rescaled variables using the tension of the D3 brane, T3, 

f = ^r, t^VTw 1 , (6) 
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we included potential terms in the brane action and we defined 

V = det (^ + / GusTdrfdytf) , (7) 

where Gu((f> K ) will play the role of a metric in the space of the scalar fields cf) 1 . By defining the mixed 
kinetic terms for the scalar fields 

X IJ = -\g^d^d^ J , (8) 

it has been shown that the explicit expression of T> reads [1] 

V = 1 - 2fG u X IJ + 4f 2 X l /xj ] - SfX [ lx J jX^ ] + 16f 4 X [ /xjX% X^ ] , (9) 

where the brackets denote antisymmetrisation on the field indices and Xj = GjkX kj . Similarly, one can 
express y/^yR[y] in terms of the fields and the geometrical quantities associated to the cosmological metric, 
leading to a multifield relativistic extension of the quartic Galileon Lagrangian in curved spacetime. 

In this paper, for simplicity, we only consider the single-field model where we can ignore the dynamics of 
angular directions and the late time curvature perturbation is dominated by the radial fluctuations. There 
are two parameters in the single filed model. One is the background value of T>: 

c 2 D ee 1 - fa 2 , (10) 



where & = y Gjjft^ plays the role of an effective collective velocity of the fields. In the DBI inflation 
model, cd corresponds to the sound speed of the perturbations, c s . The other parameter characterises the 
effect of the induced gravity 

fH 2 M 2 , . 



c 



D 



Vh 



where H is the Hubble parameter. If a = we reproduce the DBI inflation model. The parameter a is 
restricted to be < a < 1/9. 

We only show the final results for the non-Gaussianity of the gravitational potential. Detailed calculations 
can be found in Refs. |10^ 111]. The bispectrum of the Newtonian potential $ has the form 

= (2vr) 3 5^$! + k 2 + k 3 ) S(h,k 2 , k 3 ) , (12) 

where S is the primordial bispectrum shape and the Dirac delta enforces spatial homogeneity. In the 
single-field DBI Galileon model, two bispectrum shapes arise |15j : 



S^ d \h,k 2 ,k 3 ) = -^f^ d Al (13) 
(24FT 3 6 - 8K 2 2 K 3 3 Ki - SK^Ki 2 + 22K^K X Z - SK^K^ + 2Ki 6 ) 



K 3 9 Ki 
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S^ih, k 2 , k 3 ) = 162 f*£*Al • 



where 



Ki = h + k 2 + k 3 , (14) 
K 2 = (hk 2 + k 2 k 3 + k 3 h) 1/2 , 
K 3 = (hk 2 k 3 ) 1/3 , 

and we used the definition of the power spectrum 

<*(£i)$(&)> = (27r) 3 6^(k 1 +k 2 )^ . (15) 
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These two shapes arise from the two distinct cubic interactions of the comoving curvature perturbation £; 
the first shape Sfe rad ) arises from C(<9iC) 2 while the second shape S^ time ^ arises from £ 3 [T5]. The amplitudes 
of these bispectrum are determined by cd and a as 



f grad _ _^_,gmd/-, ^grad _ x 2 1 ~ a ( 9 ~ 2c D ~ 3 -^) 

/nl " 324^ V 4/ ~ l-3a(3-2cf,) ' 

f time 5 time / 1 \ timc _ 1 - 3a(5 - 2c 2 D - 4A 2 + A 4 ) 

/nl ~ 81 A \ ~ l-3a(3-2c2,) 

where A is the ratio between the angular and radial speed of sound 



(16) 



1-3.(3-24) (17) 
V l-a(5-2cf,) K ' 

3 CMB temperature bispectrum and templates 

The shapes of the bispectrum described in section 2 are not factorisable thus it is numerically challenging 
to construct optimal estimators. Instead, the WMAP collaboration has been using the following equilateral 
[T6] and orthogonal [15] templates 



S^\h,k 2 ,h) = /n1-6A|-(-^-^-^-t4t2+t4t3+ (5 perm.)), (18) 

Y rh^nj^ 1 3 2 3 1 2 3 1 2 3 / 

S^(h,k 2 M) = f^ h -^l.(-^-^-^-^+" + (5 P errn.)), (19) 

y ^ a,- 2 1 3 2 3 1 2 3 1 2 3 / 

and gave constraints on the non-linear parameters (/nl'/nl' 1 )- ^ n or der to constrain the amplitude of the 
galileon shapes in Eq. (fT3]) . we first need to relate (f^ , /n™ ) to the observed (f^, /nl 11 )- To this purpose, 
we calculate the angular bispectrum of the CMB temperature anisotropics for all the considered shapes. 
The CMB bispectrum is the 3-point function of the observed temperature anisotropies 

Wimi %m2 %m 3 ), (20) 

where a^ m are the temperature multipoles given by 

^3$(k)A,(A;)y; m (k), (21) 

and Ag is the radiative transfer function of the temperature. By using statistical isotropy, we can define an 
angle-averaged bispectrum, b^ 2 ^ 3 , as 

/ \ ( ^ £ * ^ \ R (00^ 

{ae imi a£ 2m2 ai 3m3 ) = I ^ ^ j i%^ 3 . (22) 

We compute the temperature bispectrum by projecting the primordial bispectum for the Newtonian potential 
in Eq. (JT2J) on the sky today [171 US] : 

£ Wa = Co (-|) j drr2 j dkl dk * dk * ^ k * k ^ S ( k ^ k ^ k *) ( 23 ) 
j tl (rki) A h (ki)je 3 {rk 2 ) A^ 2 (k 2 ) je 3 (rk 3 ) A t , A (k 3 ) , 
where je is the spherical Bessel function of order t and Q are the Gaunt coefficients. 



4 



The 2D scalar product of two angular bispectra, or Fisher matrix element [17], is defined as 

fl(i). B (j') = y g^Wag^Wg ; (24) 

where C\ is the angular power spectrum of temperature fluctuations, A^ 1 £ 2 £ 3 = 1,2,6 for triangles with no, 
two or three equal sides, and Z max is the maximum angular resolution attainable with the considered CMB 
survey. Note that in our analysis we include the noise contribution of the CMB survey in Cg as done, e. g., 
in Ref. [31]. 

Using the 2D scalar product, we obtain the relation between (f^ d , /n™ ) and (/nl'/nl 11 ) as P^l 

£}(grad).£j(eq) \ / g(timc) . ^j(eq) 



forth I I / £(grad).£(orth) \ / £(timc) .^g(orth) \ I I ftiiriC 

JNL / y ^ B (orth). B (orth) J y B (orth). B (orth) J ) f =1 V NL 

We should stress that these relations are obtained assuming that only one type of the templates is present 
at the same time, that is, we have either or /§L h - This is the same definition used by the WMAP team 
when they quote constraints on these non-linear parameters. 

In this paper, we are interested in a joint analysis where both equilateral and orthogonal non-Gaussianity 
exist. For the joint analysis, we introduce a new set of non-linear parameters, /* = (/nlj/nl* 1 )' which are 

■eq forth 
NL' JNL 

Y] Fijf = F u fi, i = eq, orth, (26) 



related to /< = (/° q L , f^ h ) as 



j 



where the template Fisher matrix Fi 



f B^ • B^ • B(° IthS > \ 

I B^ ■ B(° lth ^ #( orth ) • St™* 11 ) ) ^ ' 

^ ' /nl=1 

encodes the overlap between the two observational templates in ^-space. The new parameters, take into 
account the contamination from the other type of non-Gaussianity; they are equivalent to fi only if there is 
no correlation between the two estimators, that is if r = Fij/(FaFjj) 1 ^ 2 = 0. 

Using Eqs. ([25]) . ([26]) and ([2~T|). we obtain the relation between the model parameters (f^ d ,f^ e ) and 
(/nl./nl 1 ) ™M 



fcq \ / ,grad \ / n(grad) . n(eq) n(time) . n(eq) \ 

prth / _ I ftime J J I £?(grad) . _g(orth) _g(time) . _g(orth) J ' \^°) 

J NL / V J NL / \ / i /*nl = 1 

where M is the overlap matrix between the theoretical shapes and the observational templates. Given the 
best-fit values of f % from the data, f^ est , and the associated covariance matrix C = F~ , we define a x 2 
statistic for model parameters /' as 

X 2 = (f ~ fLt) (f - fLt) = (f - /best) Fij (f - fLt)- ( 2 9) 
This x 2 statistic quantifies an agreement between the observed bispectrum and the model bispectrum. 

4 Constraints from WMAP9 

In the WMAP9 paper [9], the WMAP team gave constraints on the equilateral and orthogonal non- 
Gaussianity at la level as 

/* q L = 51 ± 136, f$ L th = -245 ± 100. (30) 
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We should again bear in mind that each of these constraints are obtained assuming the absence of the other 
non-Gaussianity. The template Fisher matrix obtained from the angular bispectrum is given in the WMAP9 
paper as 



£(eq) . J3(eq) Jj(eq) . £(orth) \ / QM q 2Q 

^(eq) . ^(orth) fi (orth) . ^(orth) J ~ I q.20 1.00 



X 10~ 4 . (31) 



The WMAP9 paper also provided a relation between the theoretical parameters and (/^X' /nl 11 ) m t erms 
of the parameters of the effective field theory (c s , ^4): 



.,2 



/nl = —^(-0.276 + 0.0785A), 



1 - c 2 

/<f L th = — ^(0.0157 - 0.0163^). (32) 



In terms of (c s ,A), the theoretical amplitudes of the DBI Galileon shapes, (/pS , /nl 16 )' are gi ven by 



e gva.d _ 85 

324 V " & 



JNL on a \ 1 n 2 



iss-- -u A y-i?j' (33) 

which means that we can convert Eq. ([32]) into a relation between (f^ d ,f^ e ) and (/nL'/nl 11 ) us i n g 

(34) 



/nl A _ ( 1-052 1.271 \ / /f ad 



NL 

/ N r L th y v -°- 06 -°- 264 j v / N r 

As stressed in the WMAP9 paper, the numerical values are specific to the nine-year WMAP data as they 
are obtained using WMAP noise properties. 

We define a \ 2 statistic in terms of p = {f^, /nl 11 ) as 

X 2 = (f ~ fwMAp)^i(/ j " / W MAp)- (35) 

The central values /^map are obtained by substituting the values of Eq. (j30|) into Eq. ([26]) . and by making 
use of the template Fisher matrix in Eq. (|34p : 



/J-^j/wMAP = F H fi WMAP, fi WMAP = (51,-245) . (36) 

j 

This yields the following central values for the non-linear parameters 

(/nI /NL h ) = (153-1, -275.6) , (f^ d , f*£») = (-1537, 1393) . (37) 

In Fig. 1, the left panel shows 68%, 95% and 99.7% confidence regions in the (/nl'/nl 11 ) pl ane defined 
by threshold \ 2 values 2.28, 5.99 and 11.62. The right panel shows the same confidence regions in the 
(/nl* 1 ' /n™ 6 ) pl ane - We can convert these confidence regions to those of the theoretical parameters (co,a) 
and (c s ,A) using Eq. (fT6]) and Eq. ([33]) . respectively. In Fig 2, these constraints are shown in the (c s ,A) 
plane in the effective theory and the (cd, a) plane in the DBI Galileon model. WMAP9 data prefers non-zero 
a at 95% level because / N L h is non-zero at 95% level. However, the DBI inflation model with A = — 1 or 
a = is still allowed wihtin 99.7% level. 
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Figure 1: WMAP9 constraints on (/nl'/nl* 1 ) m the ^t P ane l an d (/nl^\ /n™ ) i n the right panel, 
(black, solid), 95% (blue, dashed) and 99.7% (red, dotted) confidence contours are shown. 



5 Constraints from Planck 



In this section, we estimate how well Planck's measurement of (/nl'^nl* 1 ) wm constrain the theoretical 
parameters in the DBI Galileon model. We utilise the Fisher module of the Second Order Non-Gaussianity 
(SONG) code [20] to numerically obtain the F and M matrices which are needed to relate (/nl'/nl 11 ) 
to /nl C ) through Eq. ([28]). In doing so, we take into account the expected sensitivity and noise 

properties of Planclfi 

In order to obtain the Fisher matrices, we first estimate the bispectrum integral in Eq. (|23p for the 
four shapes considered in this paper: equilateral, orthogonal, and the two galileon shapes in Eq. (|13|) . The 
equilateral and orthogonal templates are separable in (^1,^2,^3), meaning that their computation can be 
conveniently split into one-dimensional integrations. The galileon shapes do not have this desirable property 
and we treat them as described in Sec. 5 of Ref. [20]. We obtain the temperature transfer functions A;(/c) with 
CLASS [21] by employing a LCDM model with WMAP9 cosmological parameters [23] whereby h = 0.697, 
Q b = 0.0461, Q cdm = 0.236, n A = 0.7181, A s = 2.43 x 10~ 9 , r rcio = 0.08, N cS = 3.04. Note that, since 
the galileon shapes were computed assuming slow-roll conditions, we consistently set n s = 1 also for the 
equilateral and orthogonal templates. We checked that our results do not depend on this assumption. 

We compute the full 4x4 Fisher matrix in Eq. (|24p by interpolating our four numerical bispectra 
on a mesh in (£1,12,^3) [221 I20j . We employ the noise model described in Ref. [TO] and consider the 
100, 143, 217 Ghz frequency channels of the Planck experiment, with noise and beam parameters from 
Ref. |24j . As a result, we obtain the following full Fisher matrix 



full 



/ 2.38 -0.208 2.51 3.06 \ 

-0.208 8.47 -0.708 -2.49 

2.51 -0.708 2.68 3.39 

\ 3.06 -2.49 3.39 4.67 / 



x 10" 



(38) 



where the ordering of the rows and columns is f^, /§L h ' /nl^' /nl 16 - ^he -Ff„u matrix contains all the 



*We plan to release the code that we used for this analysis as a separate module for the Boltzmann code CLASS [21] by the 
end of 2013 
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Figure 2: WMAP9 constraints on (c s ,A) in the left panel and (a,co). Confidence contours are the same as 
in Fig. 1. 



information needed to produce Planck forecasts for the parameters of the DBI Galileon model. It also 
encodes the correlations in £-space between the two considered Galileon shapes (lower-right submatrix). 
The F matrix is the upper-left submatrix of the full Fisher matrix 

/ S (eq).£(eq) £(eq) . fl (orth) \ / 2 .38 -0.208 \ 4 

* ~ Dfenl D (orth1 oforth) oforth) ~ n one Q AT X 1U • l^ 9 ) 



fl (eq) . £(orth) £(orth) . £(orth) I _Q.208 8.47 

7 JNL = 1 v 

The errors on and /£ r L th can be obtained as l/y/Wi, giving A/^ = 64.8 and A/§ r L th = 34.4. Note that 
these constraints do not include polarisation, which will improve the bounds by a factor of 0(2) [25]. 

The M matrix, that is the overlap between theoretical shapes and templates, is the upper-right submatrix 
of the full Fisher matrix in Eq. ([38]): 

( B(sn«i) . fl(«0 B (timc) . B ( cq ) \ / 2.51 3.06 \ 4 , . 

~ V Bfe rad ) • B(° rtl1 ) #( timc ) • B(° rth ) J ~ V -0.708 -2.49 J X ' ^ ' 

^ ' /nl=1 ^ ' 

Then the relation between (/§ r L ad , e ) and (/nL/nl 11 ) is S iven b Y Ec l- dZHD as 



/n! A = ( 1-048 1-261 \ / f^ d 
Ml h J V "0-0578 -0.263 J \ 



(41) 



Although these numerical coefficients are specific to Planck, they are very similar to WMAP9 results. This 
is partly explained by the fact that the coefficients are obtained as ratios of Fisher matrix elements, i. e. they 
measure an overlap rather than an amplitude. We again define a x 2 statistic 



r = (/*-/U)W-/*J 5 ( 42 ) 

where f^ est are the best-fit values. 

We consider two hypothetical cases. In the first scenario, we assume that Planck will measure the same 
central values found in WMAP9, Eq. ([3"T|) and use them as our best-fit values. In the second case, we assume 



that Planck will find no evidence of non-Gaussianity (/nL'/nl 11 ) = (0)0)- in Fig- 3, the left panel shows 
68%, 95% and 99.7% confidence regions in the (f^ q L , f^) plane defined by threshold x 2 values 2.28, 5.99 
and 11.62. The right panel shows the same confidence regions in the (f^ d , /n™) plane. In Fig 4, these 
constraints are shown in the (c s ,A) plane in the effective theory and the (co,a) plane in the DBI Galileon 
model. In the first case, even at 99.7% level the contour closes in the (c s ,A) or (co,a) planes, significantly 
reducing the parameter space. In particular, we can rule out the DBI inflation model (A = -1 or a = 0) 
with high statistical significance (x 2 > 50). On the other hand, in the second case, we only get a lower 
bound for c s or cp. 




Figure 3: Planck constraints on (/nli/nl* 1 ) in the P ane l an(1 (/ni? > /n™ 6 ) 111 the right panel. 68% 
(black, solid), 95% (blue, dashed) and 99.7% (red, dotted) confidence contours are shown. 



6 Conclusion 

In this paper, we obtained constraints on the two parameters (a, Co) in the single-field DBI Galileon model 
using WMAP9 results. The parameter a parametrises the effect of the induced gravity on a brane and 
describes the deviation from the conventional DBI inflation model, while cd becomes the sound speed in 
the DBI inflation limit. The WMAP9 results prefer non-zero a, although the DBI inflation model is still 
allowed within 99.7% level. The bispectrum of the Newtonian potential in this model is not separable and it 
is therefore numerically challenging to construct an optimal estimator. Therefore, we used the SONG code 
[20] to obtain the relation between the amplitudes of theoretical bispectra {f^ d , /n™ ) an d the equilateral 
and orthogonal observational templates (/^l, /nl* 1 ) by properly taking into account the sensitivity and noise 
properties of Planck. We then used the bispectrum Fisher matrix for (f^, /nl 11 ) to predict the constraints 
on the two parameters. 

We found that the conventional DBI inflation model would be strongly disfavoured if the Planck ex- 
periments were to confirm the central value obtained by WMAP for (/nl'^nl* 1 )- This is driven by the 
non-zero "detection" of the orthogonal non-Gaussianity. At 95% level, we obtained the constraints as 
0.0964 < a < 0.0973 and 0.007 < cd < 0.0105. On the other hand, if the central values turn out to be 
zero, we will only obtain the lower bound for co- Note that our forecasts do not include polarisation, which 
would improve the constraints further. 
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0.100 



Figure 4: Planck constraints on (c s , A) in the left panel and (a, cd) in the right panel. Confidence contours 
are the same as in Fig. 3. 



As noticed by Ref. [26], the WMAP9 best-fit parameter for cd tends to be small because the orthogonal 
shape is generated by a cancelation of the two Galileon shapes in Eq. (|13p . This leads to a large trispectrum 
as its amplitude is proportional to l/dp. Thus the trispectrum will provide an additional consistency test 
of the model. 
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